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Abstract 

Elements of the spin density matrix for W bosons in e + e~ — > W + W~ — > qq'ivi events are measured 
from data recorded by the OPAL detector at LEP. This information is used to calculate polarised differ- 
ential cross-sections and to search for CP-violating effects. Results are presented for W bosons produced 
in e + e~ collisions with centre-of-mass energies between 183 GeV and 209 GeV. The average fraction of 
W bosons that are longitudinally polarised is found to be (23.9 ± 2.1 ± 1.1)% compared to a Standard 
Model prediction of (23.9 ± 0.1)%. All results are consistent with CP conservation. 
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1 Introduction 



We present a spin density matrix (SDM) analysis [1] of W bosons pair-produced in e + e collisions at 
LEP using data collected by the OPAL collaboration with centre-of-mass energies between 183 GeV and 
209 GeV. 

The Standard Model (SM) tree-level Feynman diagrams for W pair production at LEP are the s- 
channel diagrams with either a Z° or a photon propagator and the i-channel neutrino exchange diagram. 
Following the standard nomenclature, these three charged current diagrams are collectively referred to 
as CC03 [2]. The WWZ and WW7 vertices in the s-channel Feynman diagrams represent triple gauge 
couplings (TGC). 

By measuring the spin state of the W bosons we can investigate the physics of the TGC vertices and 
thereby test the SU(2)l x U(1)y gauge structure of the electroweak sector. The longitudinal helicity 
component of the spin is of particular interest as it arises in the SM through the electroweak symmetry- 
breaking mechanism which generates the mass of the W. In addition, comparisons of the spin states of 
the W~ and W + are sensitive to CP-violating effects. Such effects are absent both at tree-level and at 
the one-loop level in the SM. 

In order to probe the spin state of a W boson, it is necessary to reconstruct the directions in which 
its decay products are emitted. The polar angle of the W~ with respect to the electron beam direction 
is denoted by #w throughout this paper. The polar and azimuthal angles of the outgoing fermion in the 
rest frame of the parent W~ are denoted by 9^ and 4>*f respectively 1 . The polar and azimuthal angles of 
the outgoing anti-fermion in the rest frame of the parent W + are denoted by 9*j- and <p*j respectively. 

The SDM formalism is described in section 2 whilst the practical implementation is explained in 
section 4. Section 3 details the data samples and Monte Carlo simulations used in this analysis. The 
sources of systematic uncertainty in the results are explained in section 5. The measured fractions of lon- 
gitudinally polarised W bosons and other related results are presented in section 6 and their implications 
discussed in section 7. 



2 Representation of the W boson spin state 



2.1 The spin density matrix 



The angular distributions of the decay products of an ensemble of W bosons can be analysed to determine 
the composition of the average spin state, which depends on the centre-of-mass energy of the reaction, 
y/s, and the polar production angle of the W~, #w It is convenient to represent this composition by the 
spin density matrix, p TT /, defined as [1]: 



Ptt> 



(s, cos #w) = 



\,X',T 



(A,V) 



(1) 



lr The axes of the right-handed coordinate system in the parent W rest frame are defined using the helicity axes convention 
such that the z*-axis is along the boost direction, b, from the laboratory frame, and the j/*-axis is in the direction e~* x b, 
where e~* is the direction of the incoming electron beam. 
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(A A') _ 

where F T ' is the helicity amplitude for producing a W with helicity r from an electron with helicity 
A and positron with helicity A'. As the W boson has unit spin, r can take the values +1, or —1. The 
SDM for the W + is defined analogously in terms of the W + helicity amplitudes. 

The SDM is an Hermitian matrix with unit trace, and is fully described by eight free parameters. The 
elements lying on its major diagonal are the probabilities of observing a W boson in each of the three 
possible helicity states and are therefore positive as well as purely real. The real and imaginary parts of 
the off-diagonal terms measure the interference between the helicity amplitudes. 



2.2 Projection operators 

Assuming a V — A structure for the W~ coupling to fermions, the expected angular distribution for 
massless fermions in the rest frame of the parent W~ is given in terms of the diagonal elements of the 
SDM [3] by: 

3 23 
P— g (l + cosfl}) + poo- sin 2 9* f 

+/ 9 ++ jj(l-cos^) 2 . (2) 



1 dcr 

a dcos 9* f 



In this paper, a is the cross-section for the process e + e~ — ► W + W~ — > qq'£v£. In practice, the ratios of 
the masses of the SM leptons and quarks to the mass of the W boson are sufficiently small that deviations 
from equation (2) are negligible compared to the statistical precision of the measurements being made 2 . 
Hence, it is possible to extract the diagonal elements of the SDM by fitting this function to the cos dj 
distribution obtained from the data. Such an analysis has been published by the L3 collaboration [4]. 

As in previous OPAL analyses [5,6], we use the alternative method of constructing projection opera- 
tors, A TT >, which satisfy: 



f a — a d I a a*AA.* ■ A TT / dcos 9*fd6* f 

w _ J dcOS e W dcOS 0*d« TT / Yf 

Ptt> («, COS 0w) = da • ( 3 ) 

dcos 



Equations (4) to (6), listed below, give the projection operators used to extract the diagonal elements 
of the W~ SDM. They depend only on the polar production angle of the fermion, 9*j. The projection 
operators corresponding to the off-diagonal elements of the W~ SDM are given in equations (7) to (9) 
and have azimuthal angular dependence. 



1) (4) 

(5) 

1) (6) 

A+_ = 2e 2ir f (7) 

2 The top quark is too massive to be produced from on-shell W bosons and the production of bottom quarks is highly 
suppressed. 



A__ = i(5cos 2 0} + 2 cob 6} 
A 00 = 2-5 cos 2 0) 
A ++ = i(5cos 2 0*-2cos0* 
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A+o = ^-(l-4cos^)e-^ (8) 
A-o = ^-(l + 4co S ^)e^ (9) 

The operators for the W + can be obtained by replacing cos 6>j and 4>J by cos 9J and 4>J. 

2.3 Polarised cross-sections 

The longitudinally and transversely polarised differential cross-sections for W production are given in 
terms of the SDM elements by [1]: 



da L da 

Poo " 3 7T~ ( 10 ) 



dcos #w dcos 6*w 

dcT . , da 

(P++ + P-)-^znr- • ( n ) 



dcos #w dcos 6*w 

The electric charge of the charged lepton in a qq'ii' event can be reliably reconstructed so that there is 
no ambiguity in determining cos 6*j or <pj for use in the projection operators. The charge of a quark which 
originates a hadronic jet is not readily accessible, so the measured angular distributions for hadronically 
decaying W bosons are folded such that cos 8*j- lies between and 1 and 4>*f lies between and ir. Although 

neither p ++ nor p can be measured individually after this folding, their sum is unchanged and the 

polarised differential cross-sections can still be evaluated. 

The total polarised cross-sections are obtained by integrating the differential cross-sections with re- 
spect to cos 6>\y. 



2.4 Effects of CP violation 

Much of the sensitivity of W pair production to CP- violating interactions is contained in the distributions 
of the azimuthal angles <p*j and <p*j. Both azimuthal angular distributions are symmetric about zero at 
tree-level in the SM. The presence of a CP-violating phase at the TGC vertex would, in general, shift the 
distributions to introduce an asymmetry. This effect can be measured from the off-diagonal elements of 
the SDM, as described below. 

Under the assumption of CP invariance, the SDM for the W~ and the SDM for the W + are related 
by [7]: 

Ptt> — P-t-t' ■ \ LZ ) 

The time-reversal operator T can be approximated by the pseudo time-reversal operator T which 
transforms the helicity amplitudes into their complex conjugates rather than interchanging their initial 
and final states. At tree-level, the effect of the pseudo time-reversal operator is exactly equivalent to the 
effect of the true time reversal operator [8,9]. Under the assumption of CPT invariance, the SDM for the 
W~ and the SDM for the W + are related by: 
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- (^) 



(13) 



It follows that tree-level CP non-conserving effects will only violate the imaginary part of equation 
(12). This motivates the construction of a^ T , and which are defined below in terms of the imaginary 
parts of the off-diagonal SDM elements and are measured in units of cross-section: 



(7 



w- 



/r /m {^'}-d^ dcos ^ 

d<7 



-1 

H-l 



(14) 
(15) 



From these quantities we form experimentally accessible CP-odd observables: 



A™ 


= a +- 


— <7. 


a cp 

aCP 
^-0 


- rr W ~ 

— a +0 

= 0--0 


— a 

— cr. 



(16) 
(17) 

% + • (18) 



We also construct CPT-odd observables sensitive to the presence of loop effects: 

a cpt = a w- +(7 w+ (19) 

AS* = -^+- W o + (20) 

AS* = ^ 0+ aX . (21) 



No assumptions about the form of the TGC vertices are necessary to extract these observables from the 
data. Hence this study is complementary to the CP-violating TGC parameter measurements previously 
published by ALEPH and OPAL [5,10]. 



3 Data and Monte Carlo simulations 



3.1 Data sample 

This analysis used data collected by the OPAL detector [11] at LEP during the years 1997 to 2000. The 
data were collected with centre-of-mass energies clustered around eight nominal energy points: 183 GeV, 
189 GeV, 192 GeV, 196 GeV, 200 GeV, 202 GeV, 205 GeV and 207 GeV. 

Events with qq'tu final states originating from a pair of W bosons, where the charged lepton can 
be either an electron, muon or r-lepton, are referred to as signal in the remainder of the paper. The 
final state can additionally contain any number of photons. Events of all other types are referred to as 
background. Only those data events compatible with the signal definition were used to form the SDM. 
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The luminosity-weighted mean centre-of-mass energies and integrated luminosities of the data samples, 
as obtained from measurements of small angle Bhabha events in the silicon tungsten forward calorime- 
ter [12], are listed in table 1. The total integrated luminosity was 678.5 pb _1 , which corresponds to a SM 
prediction of approximately 5000 signal events being produced. The eight data samples were analysed 
separately and the results from each of these analyses were then combined using the method described 
in section 6. 



(GeV) 


jCdtipb- 1 ) 


observed events 


expected events 


182.7 


57.40 


329 


331.0 


188.6 


183.0 


1090 


1124.7 


191.6 


29.3 


170 


182.7 


195.5 


76.4 


511 


483.2 


199.5 


76.6 


457 


479.4 


201.6 


37.7 


242 


237.4 


204.9 


81.6 


482 


516.0 


206.6 


136.5 


895 


862.7 


Total 


678.5 


4176 


4214.5 



Table 1: Mean centre-of-mass energies and integrated luminosity values for the data. The number of data 
events passing the event selection detailed in section 4.1 and the expected number of events as calculated 
using the Monte Carlo simulations listed in section 3.2 are also shown. 



3.2 Monte Carlo simulations 

Samples of simulated data events with four-fermion final states consistent with having been produced 
via a pair of W bosons were generated using the KandY Monte Carlo (MC) generator formed from 
the YFSWW3 [13] and KORALW 1.51 [14] software packages. These events were weighted by factors 
calculated using KandY to provide CC03 signal samples (the CC03 Feynman diagrams are described 
in section 1). Four-fermion final states inconsistent with having been produced via a pair of W bosons 
were generated using KORALW 1.42 [15]. Additional samples of four-fermion and CC03 events used 
in studies of systematic effects were also generated using KORALW 1.42. The EXCALIBUR [16] MC 
generator was used to re- weight four-fermion events for parts of the systematic error analysis. The KandY 
and KORALW samples used in the main analysis were hadronised using JETSET [17]. To estimate the 
fragmentation and hadronisation systematic uncertainties, HERWIG 6.2 [18] and ARIADNE 4.11 [19] 
were used as alternatives to JETSET for simulating hadronisation in some KORALW samples. 

In addition to four-fermion events, only quark-pair or two-photon events were found to be significant 
sources of background for the event selection described in section 4.1. Samples of Z°/7 — > qq events were 
generated by KK2F [20] and hadronised by JETSET whilst multi-peripheral two-photon processes with 
hadronic final states (e + e~ — > e + e~77 — > e + e~qq') were simulated by HERWIG. 

The MC samples were processed by the full OPAL simulation program [21] and then reconstructed 
in the same way as the data. 
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4 Measurement of the W boson spin state 



4.1 Event selection and reconstruction 

The selection algorithm applied to the data to identify qq'lv candidate events had four parts: the pre- 
selection, the likelihood selection, the kinematic fit used to derive estimates of the momentum vectors of 
the four fermions from the W decays, and the final selection. Details of the pre-selection and likelihood 
selection have been published previously in [6,22] and were based on the 172 GeV qq'lu selection described 
in appendix A of [23]. Details of the kinematic fits and final selection are given in [24]. The number of 
events passing the full selection is shown in table 1. 

For each selected qq'£u candidate event, the fitted momentum vectors of the four fermions were used 
to calculate the production and decay angles of the W bosons required for the remainder of the analysis. 
The events were then divided into eight cos #w bins of equal width. The true number of signal events 
produced in each bin was estimated by: 




where the sum is over all data events reconstructed in the bin and p/e is a detector correction 
factor (defined in section 4.2) which allows for variations with the polar and azimuthal angles of the 
efficiency, purity and angular resolutions. The unpolarised differential cross-section needed to evaluate 
equations (10), (11), (14) and (15) was estimated by dividing these numbers by the luminosities given in 
table 1. 

Using the same notation as above, the statistical estimators for the SDM elements were given by: 

The estimators for the diagonal elements of the SDM were not explicitly constrained to lie between and 
1, and hence could have unphysical values due to statistical fluctuations in the data. In addition, the 
detector correction treated p/e as statistically independent in each angular bin. Therefore, the values of 
the estimators in different cos #\y bins were statistically uncorrelated. 

As the CP and CPT symmetry tests rely on measuring asymmetries in the azimuthal angular distribu- 
tions, the folded distributions obtained from hadronically decaying W bosons were not used in evaluating 
the off-diagonal elements of the SDM. 

4.2 Detector correction 

The MC samples listed in section 3.2 were used to estimate the efficiency and purity of the event selection 
and the angular resolution with which the directions of the W bosons and their decay products were 
reconstructed in the detector. The overall selection efficiency varied between 78% and 81%, and the 
purity varied between 93% and 95% depending on the centre-of-mass energy. In addition, the efficiency 
varied between 44% and 92%, and the purity varied between 15% and 100% in the cos0\y-cos#* plane. 
Variations of the efficiency and purity with cos 9* and Z * were smaller but significant. The angular 
resolution was dominated by the measurement uncertainties in the four-momenta of the hadronically 
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decaying W bosons, but additionally included a small contribution from the 3% of charged leptons which 
were assigned an incorrect charge (mostly in the qq'rf channel) . As the generator- level MC values of #\y 
were obtained after boosting to the centre-of-mass frame of the four fermions whereas the reconstructed 
values were measured in the lab frame, the angular resolution was sensitive to initial-state radiation (ISR) 
effects. 

After reconstruction, each data event was scaled by a detector correction factor (p/e in equation 
(23)) given by the purity divided by an efficiency- like scaling factor. The scaling factor was defined as 
the ratio of signal MC events reconstructed in a given angular bin divided by the number of signal MC 
events generated in that bin, and hence included the effects of both the efficiency and angular resolution 
under the assumption that the shapes of the MC angular distributions closely approximated those of the 
data. Approximately 16% of the MC events which passed the selection were reconstructed outside of 
the cos #w bin in which they were generated, where the exact fraction depended on the centre-of-mass 
energy. In studies using the full unfolding procedures described in [25,26], the bin-to-bin migration led 
to correlations as high as 40% between the statistical errors in neighbouring cos#w bins. As the SDM 
estimators of section 4.1 are statistically uncorrelated, special care should be exercised if using the results 
of this analysis in fitting procedures. 

The numbers of angular bins used to parameterise the detector correction for leptonically and hadroni- 
cally decaying W bosons and for diagonal and off-diagonal elements of the SDM are summarised in table 2. 
These were chosen to make best use of the MC statistics whilst reducing possible bias effects. The correc- 
tions for positively and negatively charged W bosons were combined by making use of the approximate 
CP invariance of both the SM MC and the response of the OPAL detector [27]. 



SDM elements 


Decay mode 


cos 6*w 


cos d\ 


4>\ 


cos 6* 


P , P++i Poo 




8 


20 






W^qq' 


8 






10 


P+-) P+0i P-0 




8 


5 


5 




W^qq' 











Table 2: The numbers of bins, of equal width, used to parameterise the detector correction. The hadron- 
ically decaying W bosons were not used to measure the off-diagonal SDM elements. 

4.3 Bias correction 

The method employed to compensate for detector effects, described in section 4.2, tended to bias the final 
result towards the Standard Model prediction. In order to correct for this effect in the W polarisation 
measurement, a bias correction was applied to each diagonal element of the SDM and to the cos #w 

distribution. For each cos6*\y bin, the biases in the measured values of p Q0 , p and in the number of 

reconstructed events, N, were calculated from re-weighted MC samples using a grid of 40 equally spaced 

values of p 00 and 40 equally spaced values of p distributed between and 1. The bias correction, b, 

was taken to be the average bias as shown in equation (24), where the sum is over the points on the 

grid. The vector pi represents the values of p 00 and p for the i'th re-weighted MC sample. The vector 

f represents the values measured from the data sample with error matrix R. The bias in the variable 

of interest (p 00 , p or N) is denoted by b{. The bias correction for p ++ was calculated using the bias 

corrections for p 00 and p and the normalisation constraint p ++ + p + p 00 = 1. 
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Y^ ViA Pi) -V (Pi) }) 



p(f) 
1 
2 



•6,; 



(24) 



Equation (24) was derived from Bayes' theorem assuming a uniform prior probability distribution for 
pi. The values of b for elements of the SDM typically had a magnitude less than 0.1. The uncertainty 
on each bias correction was estimated by the standard deviation of the biases, 6«. These error estimates 
were typically smaller than 0.2, and are included in the statistical errors of the results in section 6. 

As the bias varied significantly over the range of measured values spanned by the statistical errors on 
the data, the uncertainty on the correction was often larger than the correction itself. The shifts in the 
fraction of longitudinally polarised W bosons due to the introduction of the bias corrections are shown 
in table 3 where the uncertainty in each shift has been evaluated by varying the bias corrections within 
their errors. In each case, the shift due to the bias corrections is smaller than the statistical error on the 
measurement (see table 5) . Tests with MC samples were used to show that this bias correction procedure 
results in unbiased estimators which give Gaussian pull distributions. 





Shifts (%) 


(GeV) 


W -+lv 


W — > qq' 


183 


-0.8 ± 1.7 


3.0 ± 1.9 


189 


-1.1 ± 1.1 


-0.8 ± 1.2 


192 


3.7 ±2.4 


2.0 ± 1.9 


196 


-0.6 ± 1.3 


-0.5 ± 1.4 


200 


1.9 ± 1.0 


2.7 ± 2.2 


202 


-0.1 ± 1.6 


2.5 ±2.4 


205 


0.6 ± 1.5 


2.6 ±2.2 


207 


-1.0 ± 1.1 


-0.3 ± 1.6 



Table 3: Shifts in the percentage of longitudinally polarised W bosons due to the bias corrections applied 
at each centre-of-mass energy. 

The computer processing time required to extend the bias correction to the off-diagonal elements of 
the SDM was prohibitive and so the simpler method described in section 5 was used as an adequate 
alternative. 



5 Systematic errors 



There are uncertainties in the shape of the MC angular distributions due to the measurement errors on 
the parameters in the SM (such as the mass of the W boson), to the incomplete description of non- 
perturbative physics effects in MC generators and to the simplifying assumptions made in modelling 
the detector response. This leads to uncertainties in the detector corrections applied to the data and 
gives systematic errors on the final results. For each source of systematic uncertainty, the full analysis 
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of the data was repeated using a range of different MC samples (or a single MC sample re-weighted 
appropriately) to form the detector correction. The full difference between the results obtained using 
each of the detector corrections in turn was assigned as the error. The total systematic error on each 
measured quantity was calculated by summing the errors associated with each source of uncertainty in 
quadrature. The error sources are listed below and their contributions to the total systematic error for 
the luminosity-weighted average fraction of longitudinally polarised W bosons are shown in table 4. 

1. The effects of the uncertainties in the modelling of each of the different background MC samples 
were evaluated by varying the contribution from the two-photon samples by a factor of two and 
the contribution from the four-fermion background samples by a factor of 1.2. These factors were 
obtained using the method described in [22]. Samples of Z°/7 — > qq events hadronised by JETSET 
were also compared to samples hadronised by HERWIG. 

2. Unlike the MC generators used in the previous SDM analyses published by the OPAL collaboration, 
the KandY MC generator includes a full treatment of 0{a) electroweak loop corrections for the 
CC03 diagrams (using a double-pole approximation). The effect of the missing higher order correc- 
tions was estimated by re- weighting the MC samples to remove 0(a) next-to-leading electroweak 
corrections and the screened Coulomb correction. 

3. The error due to uncertainty in the modelling of the jet fragmentation was estimated by comparing 
the results obtained using WW MC samples hadronised via JETSET, HERWIG and ARIADNE. 
MC samples were generated at 189 GeV, 200 GeV and 206 GeV, and used to interpolate the 
systematic errors at all eight nominal centre-of-mass energies. 

4. The W mass obtained from Tevatron and UA2 data is 80.454 ± 0.060 GeV/c 2 [28]. The systematic 
error due to the difference between this mass and the W mass used in generating the KandY MC 
samples (80.33 GeV/c 2 ) was estimated by re-weighting the KandY samples using the EXCALIBUR 
MC generator. The LEP measurements of the W mass were not used to evaluate the error, as they 
implicitly assume that the W pairs are produced via the Standard Model mechanism. 

5. The KandY MC generator includes an 0(a 3 ) treatment of initial state radiation (ISR). The effect 
of missing higher order diagrams was estimated by re- weighting the MC samples to use an O(a) 
ISR treatment. 

6. The modelling of the OPAL detector's response to hadrons and leptons was compared to that seen 
in the data using events collected at the Z° peak. The angles and energies of the jets and charged 
leptons in the reconstructed MC were smeared to give better agreement with the Z° data. The error 
associated with the uncertainty on the smearing was evaluated by varying the smearing within its 
statistical error. A fuller account of this procedure can be found in [29]. 

7. Some areas of the W boson production and decay phase space were sparsely populated by MC and 
data events. The effect of the limited MC statistics was evaluated by smearing the efficiency and 
purity corrections by their statistical errors. 

8. Biases in the off-diagonal elements of the SDM due to the form of the detector correction were 
not explicitly removed. Instead the KandY MC samples used to calculate the detector correction 
were re- weighted to simulate anomalous values of the CP-violating TGC parameters which alter the 
shape of the <p*j and <p*j angular distributions. The TGC parameters were varied by one standard 
deviation of their measured values from the OPAL 189 GeV analysis [5]. The bias associated with 
the finite width of the angular bins into which the detector correction was divided was also taken 
into account. 
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W -► qq' 


Two-photon MC T 


0.10 


0.08 


Z U /7 qq MC r 


0.30 


0.32 


Four-fermion MC 1 


0.19 


0.03 


0(a) radiation 1 " 


0.23 


0.18 


Hadronisation"!" 


0.38 


1.23 


mw^ 


0.13 


1.07 


ISRT 


0.01 


0.01 


Detector Response 


0.19 


0.06 


MC statistics 


0.29 


0.18 


Total 


0.68 


1.68 



Table 4: The systematic errors for the luminosity-weighted average percentage of longitudinally polarised 
W bosons. The error sources appear in the same order as in the list in section 5. For the combination 
procedures described in section 6, the error sources marked by a f were considered to be 100% correlated 
between the W — > tv and W — > qq' decay modes and also among the centre-of-mass energies. All other 
sources were assumed to be completely uncorrelated. 

6 Results 

6.1 W polarisation 

Table 5 shows the fraction of longitudinally polarised W bosons measured from the data samples at 
each nominal centre-of-mass energy. The values have been corrected for the detector effects described 
in section 4.2 and the bias described in section 4.3. The values obtained at 183 GeV and 189 GeV 
differ slightly from those previously published by OPAL [5,6] due to the use of improved MC generators 
(see section 5), to the inclusion of the bias correction and to minor changes to the event reconstruction 
procedure. For comparison with the data, table 5 also gives the fraction of longitudinally polarised W 
bosons predicted by the KandY MC samples. 

The measured fractions obtained from the leptonic and hadronic decays were combined using the Best 
Linear Unbiased Estimator (BLUE) method [30] in which each source of systematic error was assumed 
to be uncorrelated with all other sources of systematic error and either 100% or 0% correlated between 
the two decay modes (see table 4 for details). The correlations between the statistical errors for the two 
decay modes were measured from the data at each centre-of-mass energy and found to have a magnitude 
of less than 10%. The systematic and statistical errors were assumed to be uncorrelated with each other. 

Following the combination of the decay modes, the results from the eight centre-of-mass energies 
were themselves combined to make a luminosity-weighted average in which the correlations between 
the systematic errors were again approximated as being either 100% or 0%. The average longitudinal 
polarisation was found to be (23.9 ± 2.1 ± 1.1)%, which is in good agreement with the KandY [13, 14] 
MC prediction of (23.9 ± 0.1)%. 

The luminosity- weighted averages of p and p 00 are shown in figure 1. The polarised differential 

cross-sections are shown in figure 2. The numerical values associated with the figures can be found in [31]. 
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Longitudinal polarisation (%) 






OPAL Data 




MC 


(GeV) 


W -+lv 


W^qq' 


Combined 


Combined 


183 


18.8 ± 10.3 ± 1.1 


31.1 ± 10.4 ±2.1 


24.8 ± 7.4 ± 1.4 


26.4 ±0.2 


189 


17.6 ±5.8 ±1.0 


21.0 ±5.8 ±1.5 


19.2 ±4.3 ±1.1 


25.6 ±0.2 


192 


56.7 ± 14.3 ±0.9 


21.7 ± 13.8 ± 1.7 


38.6 ±9.7 ±1.1 


25.2 ±0.2 


196 


16.8 ±8.3 ±0.9 


7.8 ±8.3 ±1.8 


12.4 ±5.9 ±1.1 


23.9 ±0.2 


200 


32.5 ± 8.4 ± 1.1 


35.5 ± 9.4 ± 1.6 


33.9 ±6.2 ±1.1 


23.2 ±0.2 


202 


30.4 ± 11.8 ± 1.6 


33.7 ± 12.5 ± 1.8 


31.9 ±8.9 ±1.3 


23.2 ±0.2 


205 


33.1 ±8.7 ±1.4 


35.3 ± 9.3 ± 1.9 


34.1 ±6.4 ±1.4 


22.2 ±0.2 


207 


17.5 ± 6.4 ± 1.2 


21.3 ±6.8 ±1.8 


19.2 ±4.8 ±1.3 


21.7 ±0.2 


Average 


23.6 ±2.9 ±0.7 


24.5 ± 3.0 ± 1.7 


23.9 ±2.1 ± 1.1 


23.9 ±0.1 



Table 5: The fraction of longitudinal polarisation for the leptonically and hadronically decaying W bosons 
at each nominal centre-of-mass energy after detector and bias corrections. Also shown is the combined 
result obtained from the BLUE technique [30] under the assumptions described in section 6. The values 
extracted from the data are shown with both statistical and systematic errors. The values extracted from 
the generator-level KandY MC samples are shown with statistical errors only. The last row of the table 
shows the luminosity- weighted averages. 

6.2 CP-violating effects 

The luminosity-weighted averages of the CP-odd and CPT-odd observables derived from the W~ and 
W + SDM elements in section 2.4 are shown in figure 3 as functions of cos #w The results are summarised 
in table 6. Each source of systematic error in the off-diagonal elements of the SDM was assumed to be 
100% correlated between the W~ and W + results. 

Any significant deviations from zero in the first row of table 6 would constitute an unambiguous 
signature of CP violation. Any significant deviations from zero in the second row of the table would 
show the presence of loop effects. Here, all results are consistent with the SM tree-level prediction of zero 
within the quoted errors. 





A+- ( P b) 


A+o (pb) 


A_o (pb) 


CP 


0.33 ±0.17 ±0.06 


0.09 ±0.11 ±0.04 


0.02 ±0.15 ±0.06 


CPT 


-0.10 ±0.17 ±0.06 


-0.10 ± 0.11 ± 0.04 


0.07 ± 0.15 ± 0.06 



Table 6: The luminosity-weighted average values of the CP-odd and CPT-odd observables described in 
section 2.4 measured in picobarns. The values extracted from the data are shown with both statistical 
and systematic errors. 



7 Conclusion 

We have presented a spin density matrix analysis of W bosons at LEP. The diagonal elements of the 
SDM, the differential polarised cross-sections and the fractions of longitudinal polarisation all show good 
agreement with the SM prediction for each of the analysed data samples and for their luminosity-weighted 
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average. Our result is also consistent with the fraction of longitudinal polarisation recently measured by 
the L3 collaboration at LEP using pairs of W bosons decaying to the qq'ev and qq'/j,u final states [4]. 

No evidence was found for CP violation in WW production although the statistical precision of the 
analysis was not high enough to measure loop-level effects. The absence of CP-violating effects in this 
model-independent study places loose limits on the possible values of the CP-violating TGC parameters 
(hy, Ay, g%) [3], which are constrained to be of the order A^C/cr [7]. The CP-odd observable A+H 
showed the largest deviation from zero with a luminosity-weighted average value of 0.33 ± 0.17 ± 0.06 
picobarns. Using this result we find that the CP-violating TGC parameters are expected to be less than 
or of the order ©(lO" 1 ) in the centre-of-mass energy range covered by this analysis. 
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Figure 1: Luminosity- weighted averages of the diagonal elements of the SDM as functions of cos#\y- 
The points show the data after detector and bias corrections. The inner error bars show the statistical 
uncertainties and the outer error bars show the total uncertainties including both the statistical and 
systematic contributions. The histograms show the generator-level KandY MC prediction. The figures in 

the top row show p ++ and p for leptonically decaying W bosons. The figures in the middle row show 

Poo f° r leptonically decaying W bosons and hadronically decaying W bosons separately. The figure in 
the bottom row shows p 00 with the leptonic and hadronic decay modes combined. By construction, the 
diagonal elements of the SDM are related to each other by the normalisation condition p+++p__+/>oo = 1> 
but have not been individually constrained to the physically-allowed region between zero and one. 
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Figure 2: The luminosity-weighted average polarised differential cross-sections of section 2.3, where the 
average is over the eight nominal centre-of-mass energies and over the cos #w bin width. The points show 
the data after detector and bias corrections. The inner error bars show the statistical uncertainties and the 
outer error bars show the total uncertainties including both the statistical and systematic contributions. 
The histograms show the generator-level KandY MC prediction. 
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Figure 3: The luminosity-weighted average of the CP-odd and CPT-odd observables of section 2.4, 
where the average is over the eight nominal centre-of-mass energies and over the cos 0\y bin width. The 
observables are measured in units of picobarns and shown as functions of cos #\y ■ The points show the 
data after detector correction. The inner error bars show the statistical uncertainties and the outer error 
bars show the total uncertainties including both the statistical and systematic contributions. The solid 
lines show the SM tree-level prediction. 
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